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Abstract: The system we consider here is a charged fermions gas in the 
effective mass approximation, and in grand-canonical conditions. We assume 
that the particles are confined in a three dimensional cubic box A with side 
L > 1, and subjected to a constant magnetic field of intensity B > 0. Define 
the grand canonical generalized susceptibilities Xli N > 1, as successive 
partial derivatives with respect to B of the grand canonical pressure Pl- 
Denote by P m the thermodynamic limit of Pl. Our main result is that 
Xl admit as thermodynamic limit the corresponding partial derivatives with 
respect to B of P^. In this paper we only give the main steps of the proofs, 
technical details will be given elsewhere. 
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1 Introduction and results 

In this paper, we are interested in the thermodynamic behavior of perfect 
Fermi gas in the presence of a constant magnetic field B at temperature 
T > and chemical potential /i fixed. Although the particles have electric 
charge so that they can interact with the external magnetic field, we neglect 
all self-interactions and work in the effective mass approximation. We also 
neglect the spin, since it does not change the nature of our results. 

Consider that the gas is confined in a three dimensional cubic box A of 
side L > 1, centered at the origin. The constant magnetic field is B := Be^ 
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where e 3 := (0, 0, 1) is the third vector of the canonical base of R 3 . Denote 
with c the speed of light, e the electric charge which is supposed the same for 
each particle, and define the Larmor frequency uo := (-)B > 0. We associate 
to B the magnetic vector potential a defined by: 

Ba(x):=-e 3 Ax. (1) 

The operator H(uj) := \ (— zV — uoa.) 2 is essentially self-adjoint on Cq°(1R 3 ) 
[K]. Denote by Hl{uo) the restriction of this operator to Cq°(A). The Hamil- 
tonian of our one-particle-problem is the self-adjoint extension of Hl(uo) with 
Dirichlet boundary conditions. We will use the same notation for Hi(uo) and 
for its self-adjoint extension. 

Let z = e^ 4 be the fugacity (here /3 — 1/T > 0). We will allow z to take 
complex values, i.e. z G D := C\] — oo, — e^~]. The grand canonical pressure 
Pl is then given from the grand canonical partition function 5^ by 

P L (P,z,u>) = ^\nE L (P, z,uj) = ^Tr [\n (l + ze^ H ^)] . (2) 
Define for uo > 0, 



3C 



P 00 (J3,z,w) = w(2irP) 2 ^2fs(u 



2 

fc=0 



and for uo = 0, 

P oo (/?,^,0)=/3- 1 (27r/5)-i/5( 2 ) 

where f a (z) are the standard Fermi functions (see e.g. |A-Cj ). It is proved 
in |A-Cj that Pl(/3,z,uj) admits P^fi, z,u) as thermodynamic limit in the 
following sense: for all K compact included in D, and u > 0, one has 

lim sup\P L (l3,z,u)-P oo (l3,z,u)\ = 0. (3) 

Yet it is known that if we define the grand canonical density pl(/3,z,uo) := 
(3z^-((3, z,u>) and Poo(/3, z, u) := f3z^r-{f3, z,u), one also has under same 
assumptions 

lim sup \p L (ft,z,uj) - poo(/3,2:,u/)| = 0. (4) 

The grand canonical generalized susceptibilities of a gas of fermions are de- 
fined by: 

xZ(P,z,u) z,u>), N>1. (5) 
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Notice that Xl(P, z, to) is the magnetization of the system, and xiiP, z , 
is the magnetic susceptibility. In case uj = and L —>■ oo, it is known in the 
physical literature (see [A-B-N 1] for the rigorous proof) that z,0) = 

and xto(.P, z >Q) 0- Thus at zero field, the magnetic response is quadratic. 
In case uj ^ 0, the magnetization is not zero, thus the magnetic response 
becomes linear. For N > 2, Xl would give higher corrections to the linear 
approximation. 

Our main result is that the generalized susceptibilities admit a thermo- 
dynamic limit in the following sense: 

Theorem 1.1. Fix N G N*, (3 > and uj > 0. Then for every compact K 
included in D — C\] — oo, —e^\, one has 

lim sup Wl (P,z,w) -x£(/3,2 j u>)| = 0. (6) 

A straightforward consequence of this result is that under the same as- 
sumptions: 

ViV G N*, Vm G N, lim sup 

Having uniform limits with respect to z is very useful if one wants to translate 
this type of results in the canonical ensemble. See for example [C 1] and [C 
2] for further ideas in this direction. 

Now let us mention some previous works dealing with similar problems. 
One of truly rigorous results for the case u = and N — 1,2 was given by 
Angelescu et al in [A-B-N 1]. Then Macris et al in [M-M-P] discussed the 
case when uj was arbitrary but N — 1 and \z\ < 1. In [C 1] this condition on 
z was lifted for the case of Bose statistics (but this result can be immediately 
translated for the Fermi case). Concluding, our present result settles the 
question for derivatives of all order, for all Larmor frequencies and for all 
fugacities. 

We would like to remind the reader that this paper only contains very 
basic ideas about proofs, and is mainly intended to give a detailed overview 
about the long and rather complicated technical steps that are needed. 



dz 



Xl {P,z,uj)- Xco 



dz r ' 



(P,z,uj) 



0. (7) 



2 Strategy 

Let us recall Vitali-Porter theorem IH-P 
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Theorem 2.1. Let {/l}l>i be a family of holomorphic functions on a fixed 
domain D C C. Assume that < M /or all L > 1 and all z £ D. 

Assume also the existence of a subset D' C. D having an accumulation point 
z G D, such that lim£_ i . 00 exzsfo /or eac/i z G T/ien lim^oo /^(z) 

exists everywhere in D, the convergence is uniform with respect to z in any 
compact subset of D and the limit function foo(z) is holomorphic in D. 

In our case, we have D = C\] — 00,— e^~], D' = {z G C : \z\ < 1}, 

/*(•) = x£(P,;u), and /„(•) = •,")■ 

Therefore, to prove the Theorem II .![ we need to show first the pointwise 
convergence on the unit open disk 

}imx»{l3,z,u) = x2 ) (/3,z,u), \z\ < 1, (8) 

and second, that for every compact K C D one has the uniform bound in 
L > 1 

sup|xi(/3,2,u>)| < const #,£«;). (9) 
Then Theorem 11.11 would be proven. 



3 Elements of proofs 

3.1 The pointwise limit: proof of ((HI) 

Denote by (2i(L 2 (A)), || • the Banach space of trace class operators. It 
is well known that for any u £ 1, the family of operators {Wl{/3, ou)}/3>o = 
{e~ f3HL ^ S) }fs >0 is a Gibbs semigroup [H-P] (the operators H L (cu) are self- 
adjoint, positive and {Wl{(3, ^)}/3>o C Ii(Z/ 2 (A))). On the other hand Wl 
has an integral kernel Gi(x, x', /3, uS) which is continuous on A x A with 
respect to spatial variables x and x'. The diamagnetic inequality reads as: 

1 |x-x'| 2 

\G L (*,x!,(3,u)\< e~ ~^ , x,x'GA, (10) 

which then implies that 

\\W L \\ Xl =TrW L = [ G L (x,x,/3,o;)dx< — ^ (11) 
J A (27rp)2 

Suppose z in the unit disk D', from (J2J) we have (see e.g. (2.10) and (2.11) 
in [C 1] in the case of Fermi statistics): 

PlW, z,u) = 00 J2 n - Tr ( n & w )) • ( 12 ) 

n=l 
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We denote by G^x, x', /3, u>) the integral kernel of the corresponding opera- 
tor defined on the whole space: {e~^ H ^}p>o (see (2.2) in [A-C] or (4.90) in 
[C 1]). Its diagonal is very simple and is given by 

r < r \ - 1 
oolX ' X ' P '" ;J " (27r/3)3/2 s inh(^/2)- 

We remark that this quantity is x— independent and in view of (J2J) we can 
write 

P^.z.u) = -^J2 G 00 (x,x,n/3,u;). (13) 

n=l 

We are interested in derivatives of Pl with respect to u. Due to formula 
(|12|) . these derivatives will act on the trace of the semigroup. We are thus 
motivated to study the Xi-analyticity with respect to u> of the semigroup. 
Although this result was already proven in [A-B-N 1], we state it in Lemma 
13. 1\ since we will use it later on. 

In order to do that, we need to introduce further notation. Define the 
following operators by their corresponding integral kernels: 

R 1jL (x,x.',[3,uj) := a(x) • (z'V x + ujsl(x))G l (x, x', /3, u), 
R 2 , L {x,x!,(3,u) := ia 2 (x)G L (x,x',Acu). (14) 

The operators R\,l{@,w) and R2,l(P,w) are of trace class as well as the 
operators, 

I n ,L(ii,...,i n )(J},w):= / dn / dr 2 ... \ dr n W L {(3 - n,w) (15) 

Jo Jo Jo 

~ r 2> u)R i2:L (T 2 - r 3 , ^)... J R in _ liL (r rt _i - r n , uj)R injL (r n , u) 

for n > 1 and for (ii,...,i n ) G {l,2} n . We can finally give the analyticity 
result, define 

^(»i.-.«»):=( otherwise. (16) 

Lemma 3.1. The operator-valued function R 3 uj i— > Wl{/3,uj) G X\ admits 
an entire extension to C. Fix ujq G R. For all u £ C we have 

iV=0 

= A!^(-l)" ]T ^(i!,...,^)/^^,..,^)^,^). 
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This implies in particular that the traces of the semigroup Wl which 
appear in formula ()12j) are entire functions of uj. 

Remark 3.2. It is important to notice that the expansion is not really 
convenient if one wants to prove (0). That is because the expressions 
contain at least one term as a(x) which behaves like x. So direct estimates 
show that the trace norm of d d ^N L ojq) (see e.g. [A-B-N 2J) behaves like 
L 3+N , and this is very far from the desired behavior of L 3 . It is true that when 
we look at the trace and not at the trace-norm, things are quite different. In 
[A-B-N 1] it is proved atuj = and for N = 1, 2 that due to some remarkable 
identities, the terms growing like L 3+N are identically zero. What we do next 
in our paper is to give an alternative expansion which takes care of these 
singularities for all terms at the same time. 

In order to do that, we concentrate on the kernel Gl- We remark first 
the following 

Lemma 3.3. For every uj G C, the operator Wl(/3,ou) defined by the series 
in \1T\) admits an integral kernel x', j3, uj). This kernel is defined as 

the sum of a series as in {77| ) where instead of operators we consider their 
integral kernels. Then Gl is continuous with respect to the spatial variables, 
and is an entire function of uj. In addition, for all x, x' G A fixed, one has 

dN ° L < ' r \ f9 N W L \ 

x,x,P,w) = - „ (x,x,/3,a;). (18) 



duj N y ' y V 9u N , 

Hints to the proof. Notice that d dl ^ L (x, x') is kernel's derivative, while 
d d ^i L (x, x') denotes the kernel of the trace class operator d d ^j L ■ The idea 
of the proof consists in showing that when replacing I n ^ from (JT7j) with 
its integral kernel (defined as a continuous function in x, x' G A by the 
multiple convolution (|15j)). the power series in (JT7j) converges uniformly for 
x, x' G A, and has an infinite radius of convergence. The estimates rely on 
the diamagnetic inequality (jl())l and an induction argument. □ 

Since d d ffi (x, x^ f3, uj) is continuous for (x, x') G A x A, and 9 q™n l is a 
trace class operator, its trace can be expressed as the integral of the diagonal 
of its kernel (see the remark at page 523 in [K]). We conclude that for every 
uj G C: 

qjv f fi N w \ r f) N G 

^Tr (W L ((3, uj)) = Tr (/3, uj)) = ^ (x, x, (3, uj) dx. (19) 

In the light of Remark 13. 2\ we need a different formula for the above 
kernel, so that the apparent growing terms cancel each other. This will be 
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done by using a modified perturbation theory for magnetic Gibbs semigroups. 
Previous works which dealt with similar problems are [C-N], [C 1], [B-C] and 
[N]- 

For that, we introduce the magnetic phase <fi and the magnetic flux fl 
where for x, y, z G A, and as before e 3 = (0, 0, 1), 

0(x,y) := -e 3 -(yAx), fl(x, y, z) := 0(x, y) + 0(y, z) + 0(z, x). (20) 
We have 

|fl(x,y,z)| < |x-y||y-z|. (21) 
For every n > 1 and for every points x, y 1; y n G A, we introduce 

n-l 

Fli(x,yi) = 0, Fl n (x,yi,...,y„) = fl(x, y k , y k+1 ), n>2. 

k=i 

Fix uj > 0. Consider now the bounded operators given by their integral 
kernels: 

R 1>l (x,x',(3,uj) : = a(x-x') • (iV x + ua. (x)) G L (x, x! , f3 , u) , 
R 2iL (x,x',(3,u;) := ± a 2 (x - x')G L (x, x', 0,u) (22) 

and for all x G A, n > 1 and k > 

W k L>n (x,x,(3,u;):=J2(-iy £ c](n, ...^ I dr, P dr 2 ... 

j=i (n,...^)e{i,2p Jo Jo 

p- 1 > f > /" > (j (Fl 3 (x,yi,.,y 3 ))f _ , „ , 
y rfr^y dy 1 ...j dyj J - J - G L (x, y u (3 - n, u) 

^u,L(yi,y2,ri - r 2 ,w) ...i? i ._ 1)L (y :; _i,y : ,-,T ? _i - r J -,w)i2 ij . |L (y :7 -,x,r J -,w).(23) 

By convention, in the case when k — we set 0° = 1. 

The next lemma gives a new expression for the diagonal of kernel's N-th. 
derivative with respect to uj at finite volume. 

Lemma 3.4. Fix uj > 0. Then for all x G A, and /or a// iVeN*, one has 
1 r^r* ^ 

n=l 
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Hints to the proof. The proof heavily relies on the general theory de- 
veloped in [A-B-N 2] and in [C 1], where a version of Duhamel's formula is 
written for the perturbed semigroup Wl{/3, uj) for du := uj — u>o small, here 
uq > is fixed (see Proposition 3 and formula (4.61) in [C 1]). Roughly 
speaking, one has to iterate this formula N times, and identify the term 
containing (duj) N . □ 

Based on the above formula, we can give an expression for the corre- 
sponding quantities at infinite volume: 

Lemma 3.5. Fix ujq > 0. Then for all x G A, and for all N G N* ; one has 
1 d N G OQ , . 1 d N G L 



x,x,/?,w ) = lim Ar , „ (x,x,/?,a> ) 



- - /' ; I I 

= 5^5^(-l) J Yl dr ± dr 2 ... drj dyt... 



n=l j=l he{l,2} 

\N-n 



/ «Yi 7T7 vj &oo(x, yi, /3 - n, w )-Rii,Do(yi, Y2, n - r 2 , o; ) 

■■•■Rij_i,oo(yi-i > yi > T i _i - r j ,u;o) J Ri J ,oo(yi,x,r i ,Wo) (25) 
where Rj t0 o are defined as in IHj) L = oo. 



Hints to the proof. We need to estimate (G^ — G 00 )(x, x') and (z'V x + 
cj a(x))(GL — G 00 )(x, x'). We have to take into account the walls' influence 
on the integral kernel at finite volume. We use a variant of Green's formula 
for the solutions of the heat equation inside A. Thus, we get the next result. 
Define for all x G A, 5 X : = dist(x, <9A) and M := {x G A : 5 X < 1}. Then 
one has 

| (G L - Go,) (x, x', /3, w ) |, I («V X + w a(x)) (G L - G K )(x, x', /?, w ) | 

^ exp (_^( XM(x)+XM(x , )+fflp (_|_|)) !(26) 

where c > is a constant and xm is the characteristic function of M. Then 
(j2SJ) follows after some straightforward calculations from (|TUj) . (J2EJ), PHjl . the 
estimate (4.46) in [C 1], and (J2H)- □ 

Since our aim is to prove (JBJ), we write from (|12|) and (fT3"j) 

(P L - Poo)^^^) = -—J2 [ -^ / rfx(G L -G 00 )(x,x,n/5,w).(27) 

PL 3 ^ n J A 

In order to conclude that © is true, it will be sufficient to show that 



8 



Lemma 3.6. Fix uj q > 0. Then for all N e N* , one has 

< L 2 f((3,oj ,N), (21 



f f d Gl . cr Goo R A 

J \ du N X; ^ ~ x ' ^' ^ J 



where f(-,uj ,N) is a function of ft which is polynomially bounded. 

Hints to the proof. Fix N 6 N*. One denotes by F^a^x, (3, u>o) the for- 
mula obtained by replacing in formula (|2"K|) all the spatial integrals on R 3 by 
integrals on A. To estimate the difference F iiA r(x, /3, uj ) — d Q^ N L (x, x, /3, ujq) , 
we use (EU, CDl, ISU, and the estimate (4.46) in [C 1]. "Finally we find 
that 



d N G 

dx ( F LjN (ps.,P,u ) - g ^ N L (x, x, 0, uj ] 



<L 2 h{(3^N), (29) 



where fi is polynomially bounded with respect to (3. Now we need to estimate 
the difference F LiN (~x, (3, ujq) — d d ^N c (x, x, /3, u; ). From the definition of Fl jv, 
this difference will consist with integrals as in (|25j) where in at least one of 
the spatial integrals one integrates over IR 3 \ A. Since it can also be shown 
that (|2*9*jl holds again if we replace Gl with G^, the lemma is proven up to 
the use of the triangle inequality which yields (|28j). □ 



Since we have seen in f!19|) that the derivatives with respect to uj and the 
trace commute at finite volume, formula ()27)1 and Lemma iH.fil show that for 
every \z\ < 1, the derivatives with respect to uj of (Pl — Poo) ((3, z, uj) behave 
like |j, which finishes the proof of (|S]l. 

3.2 The uniform bound: proof of flHD 

If we denote by g L ({3, r, = r, z, £) = (£ - zW L ((3, uj))~ 1 zW l {t, 
then for (3 > 0, z G if C D, uj > 0, one has (see (4.2) in [C 1] for the Bose 
case) : 

P L {(3, z,uj) = ^-J^ ln(1 e + ^ Tr 09, /?, W )) (30) 

with C C C\] — oo, — 1] surrounds the eigenvalues of zWl- In addition, one 
chooses C in order to have (£ — zWl((3,uj))~ 1 bounded for every (gC and 
z E K. We have seen that uj \— > Wl((3,uj) is Ti-analytic on C; one can also 
see that Tr(gi) is a real-analytic function of uj. In the end, the generalized 
susceptibilities are well defined as functions of z on D. We see that (j^J) will 
follow from 

1 Tr ( 9l ((3,(3,uj)) 



sup sup 

zeK £ee 



< L 3 const (/3,iT). (31) 
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Here we have a similar problem as the one pointed out for the semigroup in 
Remark 13.21 We could try to use the inequality 







N 



du N 



Tr (g L (P,P,u)) 



< 



d N 9L 



(P,P,u) 



Tl 



but the right hand side behaves like L 3+N and not like L 3 as desired. What 
we do instead is finding a Taylor expansion directly for the trace, and to give 
the right estimate for its derivatives. 

In view of developing Tr(g^) as a function of u in a small real neighbor- 
hood Q of co>o > 0, we first analyze gi in Xi(L 2 (A)). As a general rule, for 
an integral operator T(a>o) with kernel t(x, x', uq), we denote by T{uj) the 
operator which has an integral kernel given by (see also (|2()j) ) 



tfx, x', uj) :- 



t(x,x.',u ). 



Lemma 3.7. Fix ujq > and N > 1. Then for every du := to — ujq small, 
there exist N trace class operators aL, n (/3,^) = aL,n(/3, w, z, (,), 1 < n < N, 
and an operator R LjN+ i(j3,u) = Rl,n+i (Pi z, £) such that 



g L (P,P,u) 



n \ N 

' n=l 



a L>n (P,u) (32) 



+ i2 



where 



\a Ltn (p,u)\\ Tl < const(p, K,C)L 3 , u e Q, 1 < n < N 



Hi^+iOMk < |^| JV+i const(/3,K,C)L J 
Notice that the operators a^ n still depend on to. 



(33) 



Hints to the proof. We use a technique which generalizes the one developed 
in [C 1], which only worked for N = 1. Our generalization is considerably 
more involved than the original argument given in [C 1], which at its turn was 
rather lengthy. See for comparison formula (4.84) in [C 1], which corresponds 
to the case N = 1 in our lemma. Full proofs will be given elsewhere. □ 



As (I32J) is valid in Xi(L 2 (A)), we can take the trace term by term in this 
equality. This gives (see also (4.85) in [C 1]): 

Tr(g L (P,P,u)) = Tr(g L (p,P,u; )) 

N 

+ ^^"f&K^w^+Ttfe+i^w)). (34) 

n=l 
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The last technical result we need is contained in the following lemma, given 
again without proof: 

Lemma 3.8. Fix ujq > 0, N > 1 and 1 < n < N. Then there exists a family 
of u -independent coefficients {&£ n (/3, wo)}meN = {b™ n {(3, Uo, z, £)}meN; an d 
a remainder r^ l (P, oS) = r^^(P, u,z,£) such that for duj = u — cu small, 
one has 

N 

Tc(a L „(P,u>)) = ^du m bl n (P,L; )+rZ+\p,L;), (35) 

m=0 

\r^(P,u)\ < L 3 \duj\ N+1 const((3,K,C), 
\bt, n (P,u Q )\ < L 3 const((3, K,C). 

Consequently, from (|3*ljl and (fHKj) we have 

N 

09, P, co , z,0 = J2 b tn n (P> «*> z > 0, (36) 

n=l 

and (|3*T|) follows. 
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